Quasi-ballistic, nonequilibrium electron distribution in inhomogeneous semiconductor 

structures 



in 
o 
o 

(N 



OS 



Dan CsontosQ and Sergio E. Ulloa 
Department of Physics and Astronomy, and Nanoscale and Quantum 
Phenomena Institute, Ohio University, Athens, Ohio 45701, USA 

We report on a study of quasi-ballistic transport in deep submicron, inhomogeneous semiconductor struc- 
tures, focusing on the analysis of signatures found in the full nonequilibrium electron distribution. We 
perform self-consistent numerical calculations of the Poisson-Boltzmann equations for a model n+ — n~ — n+ 
GaAs structure and realistic, energy-dependent scattering. We show that, in general, the electron distribu- 
tion displays significant, temperature dependent broadening and pronounced structure in the high-velocity 
tail of the distribution. The observed characteristics have a strong spatial dependence, related to the energy- 
dependence of the scattering, and the large inhomogeneous electric field variations in these systems. We show 
that in this quasi-ballistic regime, the high-velocity tail structure is due to pure ballistic transport, whereas 
the strong broadening is due to electron scattering within the channel, and at the source(drain) interfaces. 
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As the downscaling of semiconductor devices has 
reached well into the deep submicrometer regime, the 
nonequilibrium nature of the carrier dynamics has be- 
come exceedingly important. In these systems, channel 
lengths comparable to the electron mean free path are 
possible at room temperature, and thus, ballistic and 
hot electron effects are expected to strongly influence 
the transport properties of modern devices. Although 
numerous studies of hot electron and ballistic transport 
have been performed, 1 there have been very few works 
on the microscopic carrier dynamics. It is, e.g., not clear 
how the details of the scattering at the channel inter- 
faces, as well as within the channel, affect the ballistic 
nature of the transport, what the microscopic signatures 
of quasi-ballistic transport are, and how these issues are 
responsible for limiting the current in these systems. 2,3 

A comprehensive microscopic analysis of these issues 
requires the study of the electron distribution function, 
which is a challenging task, both experimentally and 
theoretically, due to the strongly nonequilibrium nature 
of the problem. There have been a few recent theoretical 
studies of ballistic MOSFETs 2,4 in which the electron 
distribution was obtained by solution of the Boltzmann 
Transport Equation (BTE) in the collisionless limit. In 
this pure ballistic limit, according to the models, the 
electron distribution within the channel is composed of 
near-equilibrium electron distributions in the contacts, 
and thus, due to the inhomogeneous self-consistent 
field, only displays an asymmetric, pronounced high- 
velocity peak, corresponding to ballistic electrons. In 
recent experimental studies on the other hand^ the 
electron distribution in submicron, inhomogeneous 
III-V systems was measured using Raman spectroscopy, 
and was shown to display strongly broadened velocity 
distributions and interesting high-velocity tail struc- 
ture. It is clear that in realistic room-temperature 
devices, the electron scattering rate is finite and thus, 
a full solution of the BTE, including scattering is 
required. Such an analysis was in fact already initiated 
by Baranger and Wilkinspi These authors solved the 



BTE self-consistently with the Poisson equation for 
n + — n — n + GaAs structures, and found significantly 
out-of-equilibrium electron distributions in the channel 
region, and in particular, a high-velocity peak corre- 
sponding to ballistic electrons. However, the physical 
mechanisms of the strong broadening and the details of 
the electron distribution characteristics were not clearly 
resolved. In addition, the scattering introduced in the 
BTE was structureless, while the distribution function 
in inhomogeneous systems is strongly energy-dependent. 

In this paper we present a numerical analysis of 
the nonequilibrium electron distribution in submicron, 
inhomogeneous semiconductor structures. We show 
that the electron distribution function has a highly 
nonequilibrium form, is significantly broadened, and 
displays interesting structure in the high-velocity tail. 
In addition, the observed features are spatially depen- 
dent and very sensitive to the inhomogeneous electric 
field, temperature and detailed characteristics of the 
scattering, where we find that an energy-dependent 
scattering process strongly affects the details of the 
electron distribution. In addition to providing insight 
to the general quasi-ballistic transport characteristics 
in the electron distribution of submicron structures, we 
believe that our results are relevant for the recently 
observed nonequilibrium velocity distributions observed 
in the recent experiments. 

Our study is based on numerical analyses in which 
we explictly calculate the ID, steady-state electron ve- 
locity distribution function, f(x,v), through a model 
n + — n~ — n + GaAs structured We consider a nondc- 
generate system and use a numerical, self-consistent ap- 
proach based on the coupled Poisson-Boltzmann equa- 
tions, which enables us to capture essential nonequilib- 
rium and inhomogeneous transport phenomenal Scat- 
tering is treated within the relaxation-time approxima- 
tion, using realistic, energy-dependent scattering rates 
corresponding to polar-optical phonon (POP) scattering, 
which is the dominating scattering mechanism in GaAs 
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at room temperature. 



The POP scattering rate is given 
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where k = (1/eoo — 1/e), £oo is the high-frequency dielec- 
tric constant, e is the static dielectric constant, Hloq is the 
optical phonon energy and no = [exp^wo/'csT) — 
is the phonon occupation number. In Eq. the first 
term corresponds to POP absorption and the second to 
POP emission, which only occurs when e(p) > hu>o. 

In Fig. n we show the electron potential energy, elec- 
tric field and distribution function at T = 300 K, and 
V b = -0.3 V, for a 200 nm long lightly doped (10 13 cm" 3 ) 
GaAs slab sandwiched between two highly doped (10 17 
cm -3 ) long contacts. In order to highlight the effects of 
inhomogcncitics and scattering while keeping the nature 
of the scattering structureless, we first use constant scat 
tering times in our calculations (in Fig.^ r 2 = 2.5 • 10 
s), and subsequently compare these results with realistic 
energy-dependent scattering times. For comparison, the 
potential energy and electric field at Vb — are plotted in 
Fig. ^a) (dashed lines). Several direct observations can 
be made from the electrostatics depicted in Fig.^a). A 
potential barrier is formed in the n~ region due to the 
inhomogeneous doping profile that gives rise to the diffu- 
sion of electrons from the highly doped n + into the lightly 
doped n~ region. Second, as a result, a strong electric 
field on the order of 10 kV/cm is formed within the de- 
vice, even in the absence of an external field. Third, as 
a bias voltage is applied, a large portion of the voltage 
drop occurs over the submicron n~ region, giving rise to 
a strongly inhomogeneous field distribution, in contrast 
to the n + contact regions, where the field in comparison 
is very low and constant. 

In Fig. Oft)) , the calculated normalized electron distri- 
bution function, f{xi,v)vth/n{xi), where Vth is the ther- 
mal velocity, is shown at different spatial points along 
the sample for the Vb — —0.3 V case. At x\ = —0.15 /j,m, 
the electrons injected from the source display a shifted 
MB distribution. This is expected throughout the highly- 
doped contact regions, where the field is constant and low 
(« 0.1 kV/cm) and thus, the distribution function can be 
described by the linear response solution to the homoge- 
neous BTE, f(x,v) = f M B(x,v)[l - veE{x)T/k B T}. At 
X2 = —0.07 /xm, corresponding to the top of the poten- 
tial barrier, the distribution is asymmetric, showing a 
suppression of electrons in the v < tail of the distribu- 
tion. This is caused by the asymmetric potential barrier, 
which prevents left-moving electrons (v < 0) to reach xi. 
Deep in the n~ region, for x > X2, the electron distri- 
bution displays a strong spatial dependence as well as 
significant broadening. To better understand the overall 
spatial dependence of the distribution function, Fig.^c) 



shows the contour plot of f(x,v), around the interesting 
n~ region. It is clear that the electron distribution in the 
n~ region displays a distinct peak that is rapidly shifted 
toward higher velocities for increasing x > xi- This is a 
signature of ballistic transport and occurs due to the rapid 
and large (compared to fcgT) potential energy dropSA^ 
and corresponding strong inhomogeneous electric field. 
In addition to the ballistic peak, for increasing x > X2, 
the electron distribution is broadened and a low-velocity 
contribution gradually builds up until the distribution is 
dominated by a low-velocity peak, see point X5 in Fig. 
^b) and the spatial dependence in Fig. QJc). The ori- 
gin of this broadening is two-fold: The thermionically 
injected electrons at X2, with v > 0, are gradually ther- 
malized at a rate 1/r as they cross the channel, thus 
gradually "depopulating" the high-velocity peak 
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FIG. 1: (color online) (a) Electric field (thin lines) and po- 
tential energy (thick lines) profile for Vb = (dashed lines) 
and Vb = —0.3 V (solid lines), (b) normalized electron distri- 
bution calculated at the spatial points denoted by the arrows 
in (a), and (c) contour plot of f(x,v)vth/n(x) around the 



region for a n 



GaAs structure with dimen- 

17 



sions 1.5/0.2/1.5 /tm and doping profile 10 17 /10 1J /10 17 cm 
Scattering time is T2 = 2.5 ■ 10~ 13 s. 
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corresponding to ballistic electrons. Second, electrons 
backscattered at the channel-drain interface, as well as 
injected electrons with v < from a near-equilibrium 
distribution in the drain, penetrate the channel with neg- 
ative velocities, contributing to the negative and low- 
velocity end of the electron distribution. These electrons 
are scattered at a rate 1/r, but also experience scattering 
at the steep potential energy barrier close to the source. 
Hence, close to the source and deep in the n~ region, 
the distribution is dominated by the ballistic peak cor- 
responding to the source-injected ballistic electrons [see 
Fig. Bib)]. Close to the drain, however, where the ef- 
fective potential barrier is lower, the drain-injected and 
reflected electrons dominate the distribution and obscure 
the ballistic peak. 

The spatial dependence of the electron velocity distri- 
bution is strongly influenced by the characteristics of the 
external and built-in fields, mobility, scattering mecha- 
nisms and temperature. In Fig. Ufa), (b) we show the 
distribution function calculated for the constant scatter- 
ing times n = 1.0 • 10~ 13 s and r 3 = 4.5 • 10~ 13 s, for 
the same structure and system parameters studied in Fig. 
^ (with T2 = 2.5 • 10~ 13 s). While energy-independent, 
these scattering times correpond to realistic mobilities of 
GaAs. For the short scattering time, n = 1.0 • 10~ 13 
s, the ballistic peak structure seen in Fig. ^c) is signif- 
icantly less pronounced and the average velocity of the 
electrons is strongly reduced. The latter observation is 
a natural consequence of the increase of the scattering 
rate. However, in addition, the fraction of the potential 
drop that occurs over the channel region also decreases 
with an increase in the scattering rate (not shown here) 
and hence, due to the decreased effective potential bar- 
rier, the drain-injected and reflected electrons penetrate 
deep into the channel. Therefore, the resulting distribu- 
tion, although strongly out-of-equilibrium, does not dis- 
play distinct ballistic features, but becomes broadened, 
asymmetric and shifted. For T3 = 4.5 • 10~ 13 s [Fig. EJb)] . 
the opposite effects are observed: A very distinct ballis- 
tic peak dominates the distribution function deep into 
the n~ region close to the drain, and less broadening is 
observed. This is due to an increase in the electron veloc- 
ity given by the decreased scattering rate, and an increase 
in the effective potential barrier seen by electrons with 
v < coming from the drain. 

In realistic systems, the scattering and relaxation pro- 
cesses can have a strong energy-dependence and thus, it is 
interesting to consider the effects of an energy-dependent 
scattering rate on the nonequilibrium distribution func- 
tion. In Fig. |2c) the energy-dependent POP scatter- 
ing time T p (e), calculated from Eq. at T = 300 K, 
is shown. For clarity, in the figure, we have also indi- 
cated the values of the constant scattering times used in 
the calculations shown in Figs.|2[a), (b) and^c). Be- 
low the threshold for POP emission, the scattering time 
is in the range « 3.5 — 4 ■ 10 -13 s, whereas above the 
threshold, the scattering time has a relatively weak de- 
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FIG. 2: (color online) Normalized distribution function, 
f{x,v)v th /n(x), for n = 1.0 ■ 10~ 13 s (a), r 3 = 4.5 ■ 10~ 13 
s (b), r p (T = 300 K) (d), and t p (T = 200 K) (f). (c) POP 
scattering time at T = 300 K. (e) Difference between the nor- 
malized distribution functions calculated for t p (T = 300 K) 
and T3 [plots (d) and (b)]. 



pendence on energy (in this particular energy interval) 
and is close to « 1 • 10" 13 s. The corresponding elec- 
tron distribution function, calculated at T = 300 K is 
shown in Fig. 0[d) . Strong signatures of quasi-ballistic 
transport are seen, i.e., formation of a ballistic peak that 
evolves toward higher velocities, similarly to the calcula- 
tions performed with the constant, energy- independent 
scattering times, and significant broadening and low- 
velocity contributions due to the scattering processes dis- 
cussed above. The details of the distribution function 
and its spatial dependence are, however, different. In 
Fig-EIe), the difference between the normalized functions 
f[x,T p (e)]v th /n(x)-f(x ) T 3 )v th /n'(x) [n(x) and?i'(:z;) are 
the electron densities corresponding to the two calculated 
distribution functions] is shown. The observed difference 
can be explained as follows: Electrons with velocities be- 
low v tp ~ 4.3 • 10 cm/s, which corresponds to a kinetic 
energy comparable to the threshold for POP emission 
[Fig.[2c)], have a large mean free path and behave ballis- 
tically, as concluded by the results shown by the calcula- 
tions for T3 shown in Fig. Gib) . For electrons with higher 
kinetic energies on the other hand the scattering rate is 
higher and thus the electrons relax toward lower ener- 
gies at an increasing rate. However, the picture is more 
complicated, since the detailed energy-dependence of the 
scattering effectively also changes the characteristics of 
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the inhomogeneous field due to the charge redistribution 
caused by the scattering. 

Finally, we briefly discuss the effects of temperature, 
which is the parameter that, next to the applied electric 
field, is most easily tuned in experiments. In Fig. [5JT) 
we show the electron distribution calculated at T = 200 
K. Below the threshold for POP emission, the calcu- 
lated scattering time (not shown here) is in the range 
rs 8.2 — 9.2-10~ 13 s, whereas above threshold, the scatter- 
ing time has weak energy dependence and is ss 1.5 • 10~ 13 
s. As anticipated from the above discussion, the main 
effects of lowering the temperature and, thus, increasing 
of the phonon scattering time, are: i) higher average ve- 
locities, ii) larger spatial extension of the ballistic peak, 
Hi) narrowing of the ballistic peak in the n _ region as 
well as the diffusive peaks in the source/drain, and iv) 
suppression of the drain-injected electron contribution to 
the electron distribution in the channel region, due to an 
increase in the effective potential barrier. 

In summary, we have studied the electron distribu- 
tion in inhomogeneous, deep submicron semiconductor 
structures by self-consistent calculations of the semiclas- 
sical BTE. We have shown that the electron distribution 



in general is strongly out-of-equilibrium, significantly 
broadened and displays pronounced structure in the 
high-velocity tail. These characteristics of quasi-ballistic 
transport are very sensitive to the energy-dependent scat- 
tering in the channel and at the source(drain) interface, 
as well as the strongly inhomogeneous electric field. Our 
results are similar to recent experimental observations. 

Note added after submission: An interesting study 
of quasiballistic transport in nanoscale semiconductor 
structures with focus on the scattering and the mathe- 
matical nature of the BTE at the top of the potential 
energy barrier at the source-channel interface has been 
published very recently by Sanopii^ We note that 
although the calculations were done for completely 
different material (Si) and system parameters, the 
features of the distribution shown in Fig. 3 of Ref. 
are similar to the features of the distributions calculated 
here, and are due to the mechanisms discussed in our 
paper. 

This work was supported by the Indiana 21st Century 
Research and Technology Fund. 
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